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Power  production  was  traditionally  dominated  by  monopolies.  After  a  long  period  of  research  and 
organisational  advances  in  international  level,  electricity  markets  have  been  deregulated  allowing 
customers  to  choose  their  provider  and  new  producers  to  compete  the  former  Public  Power  Companies. 
Vast  changes  have  been  made  in  the  European  legal  framework  but  still,  the  experience  gathered  is  not 
sufficient  to  derive  safe  conclusions  regarding  the  efficiency  and  reliability  of  deregulation.  Furthermore, 
emissions’  trading  progressively  becomes  a  reality  in  many  respects,  compliance  with  Kyoto  protocol’s 
targets  is  a  necessity,  and  stability  of  the  national  grid’s  operation  is  a  constraint  of  vital  importance. 
Consequently,  the  production  of  electricity  should  not  rely  solely  in  conventional  energy  sources  neither 
in  renewable  ones  but  on  a  mixed  structure.  Finding  this  optimal  mix  is  the  primary  objective  of  the 
study.  A  computational  tool  has  been  created,  that  simulates  and  optimises  the  future  electricity 
generation  structure  based  on  existing  as  well  as  on  emerging  technologies.  The  results  focus  on  the 
Greek  Power  Sector  and  indicate  a  gradual  decreasing  of  anticipated  C02  emissions  while  the  socio¬ 
economic  constraints  and  reliability  requirements  of  the  system  are  met.  Policy  interventions  are 
pointed  out  based  on  the  numerical  results  of  the  model. 
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Nomenclature 

* 

the  decimal  point  indicator  in  numerical  formats 

• 

the  1000  separator  in  numerical  formats 

a 

elasticity  formula  base 

&l,i 

capacity  factor  for  each  plant  type  (i) 

&m,i 

availability  factor  for  each  plant  type  (i) 

b(i ) 

learning  rate  of  plant  (i)  construction 

Q> 

capacity  orders  realised  in  the  past  (year  v  before 
the  initiation  of  the  investigated  time  period)  for 
plants  (i)  (MW) 

Cfomiv 

fixed  operational  and  maintenance  cost  of  electric¬ 
ity  production  for  plant  type  (i)  for  orders  made  in 
year  (v)  (€/MW) 

CfueliyZ 

unit  fuel  cost  of  electricity  production  during  year 
(z)  for  plant  type  (i)  ordered  in  year  (v)  (€/MWhei) 

C/,> 

investment  cost  for  orders  of  plant  type  (i)  made  in 
year  (v)  (€/MW) 

CONV 

number  of  electricity  production  technologies 
based  on  conventional  energy  sources 

dWt 

Wiener  (Brownian  motion)  vector  differential, 
normally  distributed:  ~gN  (0,1) 

dyn, 

the  natural  resource’s  maximum  potential  for  year 

(i)  (MW) 

dz 

total  demand  (including  peak  loads)  of  electricity 
for  year  (z)  (MWh) 

E 

elasticity  constant 

efft 

efficiency  of  electricity  production  for  plant  type  (i) 

(%) 

EmCo2i 

emissions  of  C02  for  a  specific  technology  (tn) 

6Tnmix 

emissions  factor  from  the  current  conventional 
generating  mix  (tn  C02/MWh) 

fco  2, 

C02  factor  of  fuel  type  (i)  (tn  C02/MWh  fuel) 

I 

total  number  of  studied  electricity  production 
technologies 

ir 

inflation  rate  (%) 

leadi 

time  required  for  the  construction  of  a  power  plant 
type  (i)  (years) 

Li,z,v 

installed  capacity  during  year  (z)  for  plant  type  (i) 
ordered  in  year  (v)  (%)  (MW) 

life, 

economic  life  for  each  plant  type  (i)  (years) 

MAD 

mean  absolute  deviation  of  the  forecast 

N  (0,1 ) 

range  of  values,  normally  distributed 

Pco2 

price  of  C02  allowances  (€/ tn  C02) 

PPz 

peak-power  demand  for  year  (z)  (MW) 

Pz 

electricity  selling  price  to  the  grid,  for  year  (z)  (€/ 
MWh) 

r 

interest  rate  (%) 

REN 

number  of  electricity  production  technologies 
based  on  renewable  energy  sources 

rm 

generation  reserves  margin  factor 

XifV  orders  made  in  year  (v)  for  each  plant  type  of 
energy  technology  (i)  (MW) 

Yr  the  entire  time  period  required  for  the  calculations 

(years) 

g  a  very  small  number  approaching  infinitesimal 

values 

6itZtV  annual  usage  factor  during  year  (z)  for  plant  type  (i) 
ordered  in  year  (v) 
ti  mean  value 

a  standard  deviation 


1.  Introduction 

The  electricity  sectors  of  many  countries  have  faced  numerous 
changes  in  their  structure  and  their  business  environment  during 
the  last  years.  First  of  all,  the  electricity  markets  have  gone  through 
a  deregulation  process,  which  has  introduced  competition  in  a 
formerly  state-regulated  economic  sector.  Therefore,  the  planning 
for  new  power  plant  additions  and  existing  plant  replacements  has 
shifted  its  focus  from  strategic  fuel  selection  to  economic 
considerations,  such  as  the  minimisation  of  the  production  cost. 
This  shift  has  also  amplified  the  effects  of  uncertainty  in  fuel  prices, 
since  now  its  effects  are  even  more  crucial  for  an  investment 
decision.  The  result  of  these  changes  is  that  the  traditional  ways  of 
deterministic  financial  analysis  of  investments,  like  the  NPV  (Net 
Present  Value)  method,  are  not  capable  of  handling  adequately  the 
increased  uncertainty  [1,2]. 

Another  significant  change  has  been  the  Kyoto  protocol. 
According  to  this  protocol,  all  the  developed  countries  that 
participate  have  committed  to  reduce  their  greenhouse  gas 
emissions  to  a  certain  target  level.  One  of  the  mechanisms 
allowing  individual  industries  to  meet  their  targets  is  the  creation 
of  emissions  trading  markets,  where  the  owner  of  emission 
allowances  may  trade  them  at  the  current  price  that  is  settled  by 
the  laws  of  demand  and  supply,  like  other  commodities.  This 
mechanism  is  of  high  importance  for  renewable  energy  projects,  as 
it  may  constitute  a  new  revenue  stream  that  will  improve  their 
financial  yield  and,  therefore,  attractiveness  [3].  On  the  other  hand, 
electricity  generation  using  fossil  fuels  is  a  major  greenhouse  gas 
emitter  and  was  one  of  the  first  business  sectors  to  receive 
restrictions  to  its  emissions.  Expectations  about  future  greenhouse 
gas  allowance  prices  already  influence  current  decision  making, 
especially  in  the  energy  sector,  which  is  greatly  affected.  As  the 
future  might  bring  a  requirement  for  purchasing  allowances  for  all 
the  greenhouse  gas  emissions  of  the  fossil-fueled  electricity 
generation,  the  uncertainty  and  volatility  of  the  greenhouse  gas 
allowances  will  play  a  critical  role  in  the  investment  decision  for 
new  power  plants. 

The  present  study  attempts  to  organise  the  above  mentioned 
considerations  in  a  computational  model  whose  objectives  are:  (a) 
to  simulate,  correlate  and  forecast  the  variables  under  uncertainty 
including  fuel  prices,  electricity  demand  and  prices  as  well  as  the 
greenhouse  gas  allowances,  (b)  to  optimise  the  structure  of 
electricity  energy  production  taking  into  account  the  stability  of 
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the  grid  as  well  as  the  Kyoto  protocol’s  conditions  and  the 
accompanied  EU  directives,  (c)  to  analyse  the  cost  and  revenue 
structure  of  the  optimised  power  sector  and  identify  the  part  due 
to  greenhouse  gas  allowances,  and  (d)  to  investigate  the  share  of 
electricity  production  based  on  renewable  energy  sources  in  the 
optimised  power  sector  and  calculate  the  expected  C02  emissions 
as  a  function  of  time.  Among  the  milestones  of  the  study  was  the 
stochastic  correlation  of  the  variables  as  well  as  the  calculation  of 
the  optimal  usage  of  the  power  plants  contributing  to  electricity 
production.  In  order  to  accomplish  the  above  milestones  and 
objectives,  some  computational  techniques  had  to  be  implemen¬ 
ted,  to  overcome  the  difficulties  arising  from  nonlinear  interaction 
of  the  various  stochastic  variables.  The  model  has  been  imple¬ 
mented  for  the  case  study  of  the  Greek  Power  Sector.  However, 
other  case  studies  are  possible,  provided  that  historical  data  of  fuel 
and  electricity  prices  are  available  and  adjusted  to  those  particular 
cases. 

The  study  is  structured  in  brief  as  follows:  in  Section  2,  a 
literature  review  is  presented.  In  Section  3  an  overview  of  the 
model  is  given  while  in  Sections  4  and  5,  details  of  the 
mathematical  and  computational  model  are  provided.  Section  6 
includes  the  results  of  the  model’s  implementation.  Finally,  in 
Section  7  a  critical  analysis  is  presented  with  suggestions  on 
strategic  interventions  and  policy  issues  related  with  deregulated 
electricity  production. 

2.  Related  studies 

In  the  field  of  power  generation  investments,  even  recently, 
traditional  cash-flow  methods  of  analysis  were  used  to  address 
power  sectors’  structure  [4].  However,  stochastic  analysis  experi¬ 
ence  increasing  acceptance  among  the  energy-related  scientific 
community.  The  answer  to  the  question  “invest  or  wait”  is  dealt  in 
recent  studies  as  well  as  in  the  present  study.  One  may 
discriminate  two  separate  classes  of  models:  (i)  those  optimising 
policy  interventions  in  future  generation  mix  thus  suggesting 
state-originating  power  licensing  procedures  [2,5]  and  (ii)  those 
who  suggest  investment  options  for  private  energy-related 
businesses  [6-8].  The  present  study  should  likely  be  classified  to 
the  first  (i)  category  of  modelling. 

Robust  models  are  developed  in  the  above  mentioned  studies  of 
category  i  (optimised  portfolios)  although  an  optimal  usage  of 
power  plants  might  be  also  incorporated,  i.e.  not  only  finding  the 
required  annual  capacity  investments  but  also  calculating  their 
optimal  capacity  factor  depending  on  the  generation  mix  over 
time.  Recently,  the  uncertainties  arising  from  climate  changes 
seem  to  gain  increasing  focus  [9,10].  These  studies  address  the 
major  issues  of  emissions  trading  that  emerge  in  parallel  with 
severe  climate  changes.  Despite  of  the  excellent  analyses,  either 
the  C02  trading  costs-revenues  are  not  calculated  inside  the 
optimisation’s  objective  function  or  there  is  no  power-portfolio’s 
optimisation  included.  In  Refs.  [11-13]  methodological,  organisa¬ 
tional  and  financial  issues  addressing  the  transition  phase  to 
deregulated  electricity  markets  can  be  seen,  but  an  analysis  based 
on  the  environmental  issues  as  well  as  emissions’  trading  issues,  is 
not  included.  Concerning  power-portfolio  optimisation  [5,14] 
integrated  multi-objective  models  have  been  developed.  An 
elaborate  consideration  of  heat  recovery  in  different  power 
generation  methods  is  included,  thus  addressing  the  combined 
heat  and  power  (CHP)  industry  too.  Cost  killing  policies  in  different 
geographical  regions  are  also  addressed.  In  Ref.  [15]  discrete 
probabilistic  optimisation  models  have  been  analysed  for  the 
evolution  of  stochastic  variables.  The  models  mentioned  above, 
were  treated  by  either  dynamic  or  stochastic  programming 
techniques.  The  SQP  (Sequential  Quadratic  Programming)  ap¬ 
proach  of  the  present  study  is  a  consequence  of  the  built-in 


interaction  of  the  various  contributing  stochastic  variables  which 
in  turn  implied  nonlinear  considerations  in  the  optimisation 
algorithm. 

3.  Methodology 

3.2.  An  overview  of  the  model 

The  electricity  market  is  assumed  to  be  deregulated,  as  far  as 
the  electricity  prices  are  concerned.  Nonetheless,  the  annual  upper 
bounds  of  licence  provisions  from  each  competing  technology  are 
assumed  to  be  controlled  by  the  state,  thus  reflecting  the  current 
situation  in  the  Greek  Power  Sector.  A  random-walk  procedure  is 
used  in  order  to  forecast  the  future  electricity  demand  and  forward 
prices.  The  same  stands  for  fuel  prices  as  well  as  for  C02  allowance 
prices.  Appropriate  stochastic  differential  equations  (SDEs)  (either 
Geometric  Brownian  Motion  (GBM— [16])  or  Mean  Reverting  (MR- 
117]))  are  numerically  solved  using  an  Euler-Marujama  solver.  The 
uncertainty  of  the  algorithm  is  considered  by  utilising  a  Monte- 
Carlo  iterative  sub-routine  embedding  that  solver.  The  next  step  of 
the  study  is  to  reconstitute  the  future,  in  an  optimised  form:  the 
previously  forecasted  data,  feed  an  optimisation  computational 
procedure  which  reproduces  the  optimal  future  structure  of  the 
power  sector,  i.e.  what  capacities  from  each  type  should  be  ordered 
yearly  and  what  should  their  cumulative  energy  production  be  for 
the  next  years.  Optimality  is  determined  when  the  maximum  Net 
Present  Value  (NPV)  of  the  system  is  achieved  while  simulta¬ 
neously  meeting  the  total  power  demand.  It  is  important  to  note 
that  the  “system”  mentioned  above  includes  the  state,  the  grid,  the 
society  and  every  power  producer  contributing  to  electricity 
production.  The  state  in  turn  includes  the  regularity  authorities 
and  the  system  operator  which  are  presented  in  Appendix  A.  The 
optimisation  procedure  comprises  of  a  nonlinear  algorithm,  based 
on  the  solution  of  the  Karush-Kuhn-Tucker  (KKT)  equations  which 
are  approximated  with  Sequential  Quadratic  Programming  (SQP). 
Investment  costs  are  calculated  using  learning  curves  of  expertise- 
cost  relationship  in  a  similar  vein  to  [2].  Emission’s  trading  costs 
and  related  revenues  are  embedded  in  the  optimisation  procedure 
as  contributors  to  the  objective  function.  Finally,  the  Kyoto 
protocol  requirements  and  the  related  EU  directives  are  also 
considered  as  basic  problem  constraints. 

3.2.  Stochastic  simulation 

The  simulation  of  stochastic  variables  is  required  for  represent¬ 
ing  their  evolution  over  time.  Several  options  are  potentially 
available:  the  class  of  ARMA  (Auto  Regressive  Moving  Average) 
models  is  mainly  useful  for  stationary  time-series  simulations  and 
the  ARIMA  (Auto  Regressive  Integrated  Moving  Average)  models 
which  are  used  for  non-stationary  time-series  [18].  The  GARCH 
(generalised  autoregressive  conditional  heteroskedasticity)  type  of 
models  [19,20]  on  the  other  hand  is  mainly  indicated  for  time- 
series  with  extremely  high  volatility  levels  (noise). 

In  the  present  article,  futures  contracts  (i.e.  using  forward 
prices)  are  rather  simulated  instead  of  spot  prices  thus  represent¬ 
ing  average  prices  over  a  time  period  which  do  not  contain 
information  about  the  very  short  term  variations  of  electricity 
prices  (as  stated  in  Ref.  [21  ]).  By  utilising  a  probability  distribution 
fitting  tool,  it  has  been  validated  that  the  time-averaged  forward 
prices  as  well  as  the  other  stochastic  variables  samples  are 
normally  distributed. 

Assumption  1.  Since  normally  distributed  samples  are  required 
for  GBM  or  MR  simulations,  it  may  be  assumed  that  the  time- 
averaged  stochastic  variables  follow  a  random-walk  process:  (i) 
the  Ito  process  (GBM  stochastic  differential  equation),  when  the 
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variables’  history  fail  to  define  a  long  run  mean  value  (electricity 
demand  and  electricity  and  fuel  prices)  and  (ii)  the  Ornstein- 
Uhlenbeck  (MR  stochastic  differential  equation)  in  the  opposite 
case  (C02  allowance  prices).  In  the  latter  case,  the  long  run  mean 
represents  market  equilibrium.  Therefore,  in  the  present  article 
forward  prices  are  modelled  with  a  GBM  process.  Both  processes 
(GBM/MR)  may  be  considered  as  particular  cases  of  ARMA  (or 
ARIMA)  models  when  dealing  with  discrete  time  and  they  are 
mathematically  formulated  with  Eqs.  (2*)  and  (3*),  respectively, 
presented  in  Appendix  B. 

Concerning  forward  electricity  prices,  Keppo  and  Lu  [6]  also 
adopted  the  above  assumption  (GBM  representation  of  forward 
price  evolution).  In  addition,  they  assumed  that  jump-processes 
are  absorbed  by  averaging  spot  prices  over  time  and  using  the 
resulting  forward  prices  in  futures  contracts.  Barlow  [22] 
developed  GBM  and  MR  derivatives  for  the  evolution  of  forward 
prices  and  their  correlation  with  spot  prices,  accounting  the 
spike  absorption  effects.  Audet  et  al.  [23]  Laurikka  [7]  Laurikka 
and  Koljonen  [8]  and  Kumbaroglu  et  al.  [2],  have  also  adopted 
various  types  of  GBM  or  MR  models  for  the  projection  of  forward 
and  spot  electricity  prices  as  well  as  for  the  evolution  of  other 
stochastic  commodities  like  fuel  prices  and  electricity  demand 
loads. 

In  the  present  work,  the  Brownian  differentials  of  the  stochastic 
variables  were  correlated  based  on  past  data  in  an  attempt  to 
capture  endogenous  dynamics.  A  further  analysis  on  the  endoge¬ 
nous  interaction  of  electricity  prices  with  other  stochastic 
variables  has  not  been  attempted.  It  would  rather  raise  challenges 
in  terms  of  computing  time  and  input  fitting,  considering  that  a 
multiple  investors  problem  is  dealt  with  many  individual  power 
plants  belonging  to  10  discrete  power  generation  technologies. 
Moreover,  sufficiently  detailed  data  would  be  needed  reflecting  the 
dynamics  of  price  variations  in  relationship  with  capacity 
instalments,  start-ups  or  shut  downs  for  all  the  participating 
technologies.  However  such  detailed  data  were  not  available. 
Nonetheless,  the  definite  impact  of  an  individual  power  investor’s 
actions  on  electricity  prices  may  be  minimised  when  multiple 
investors  are  considered  and  also,  provided  that  the  following  two 
assumptions  are  additionally  made: 

Assumption  2.  Market  equilibrium  is  considered 

Assumption  3.  The  investors’  profit  is  a  direct  function  of  electric¬ 
ity  price 

As  documented  in  Ref.  [21],  competitors’  actions  may  affect  an 
investor’s  decisions  in  two  ways:  (i)  the  option  value  of  waiting 
may  be  reduced  as  this  may  force  other  competitors  to  enter  the 
market  prior  to  the  first  investor  thus  capturing  its  potential  share 
on  the  market  and  (ii)  the  upper  limit  of  profit  distributions  may  be 
reduced  since  new  plant  instalments  usually  contribute  to  price 
dropping.  Which  effect  is  stronger  depends  on  the  characteristics 
of  price  modelling  and  on  the  other  investors’  strategies.  For 
multiple  investor  models  -  in  which  Assumptions  1,  2  and  3  are 
valid  -  the  two  effects  are  cancelled  out  and  therefore  the  impact 
on  electricity  prices  is  minimised  while  the  investment  strategy 
remains  the  same  as  if  other  investors  were  not  represented  in  the 
model  (as  proved  by  Dixit  and  Pindyck  [24]).  Thus,  in  the  present 
work,  being  a  multiple  investors  model,  the  endogenous  modelling 
of  electricity  prices  is  solely  based  on  the  Assumptions  1,  2  and  3 
and  on  their  statistical  correlation  with  the  other  participating 
stochastic  variables.  They  are  all  simulated  through  a  Monte-Carlo 
process  [16,25]  which  is  comprised  of  two  steps:  (a)  producing 
multiple  SDE  solutions  using  an  Euler-Marujama  solver  and  (b) 
averaging  the  multiple  solution  paths  in  order  to  supply  the 
optimisation  algorithm  with  the  required  input  data. 


3.3.  Optimisation 

The  general  scope  of  SQP  numerical  algorithm  is  to  approximate 
an  objective  function  with  an  easier  one.  This  sub-problem  can  then 
be  resolved  and  used  as  the  basis  of  an  iterative  process.  A 
characteristic  of  a  large  class  of  older  methods  was  the  transforma¬ 
tion  of  the  constrained  problem  to  a  basic  unconstrained  problem  by 
using  a  penalty  function  for  constraints  that  are  near  or  beyond  the 
constraint  boundary.  The  constrained  problem  is  then  solved  using  a 
sequence  of  parameterised  unconstrained  optimisations,  which 
converge  to  the  constrained  problem.  These  methods  are  now 
considered  relatively  inefficient  and  have  been  replaced  by  methods 
that  have  focused  on  the  solution  of  the  Karush-Kuhn-Tucker  (KKT) 
equations.  The  KKT  equations  are  necessary  conditions  for 
optimality  for  a  constrained  optimisation  problem.  If  the  problem 
is  a  so-called  convex  programming  problem,  meaning  that  it  has 
convex  objective  and  constraint  functions,  then  the  KKT  equations 
are  both  necessary  and  sufficient  for  a  global  solution  point. 
Concerning  the  numerical  solution  of  KKT  equations,  SQP  methods 
are  mostly  utilised  in  nonlinear  programming.  Schittkowski  [26]  has 
implemented  and  tested  an  SQP  version  that  outperforms  every 
other  tested  method  in  terms  of  efficiency,  accuracy,  and  percentage 
of  successful  solutions,  over  a  large  number  of  test  problems.  Based 
on  previous  works  [27-29],  the  method  calculates  the  Lagrangian 
factors  of  the  KKT  equations.  Each  iteration  includes  an  approxima¬ 
tion  of  the  Hessian  matrix  of  the  Lagrangian  function  using  a  quasi- 
Newton  updating  method.  This  is  then  used  to  generate  a  QP  sub¬ 
problem  whose  solution  is  used  to  form  a  search  direction  for  a  line 
search  procedure.  An  overview  of  the  SQP  method  is  found  in  [30- 
33].  During  the  last  years,  great  advances  have  been  made 
concerning  the  realisation  and  implementation  of  the  above 
mentioned  works.  More  recently,  together  with  the  development 
of  highly  efficient  processors  (CPUs),  appropriate  software  have 
been  available  which  are  able  to  reach  globally  optimum  solutions 
while  maintaining  high  standards  of  accuracy  as  well  as  computa¬ 
tional  efficiency  [34]. 

3.4.  Contribution  of  the  study 

Apart  from  the  results  produced,  which  focused  on  the  Greek 
Power  Sector,  the  methodological  contribution  of  the  study  may  be 
summarised  to  the  following  combination  of  key-points: 

(1)  Optimisation  of  the  power  sector’s  future  structure  using  a 
nonlinear  SQP  approach  able  to  improve  -  on  average  -  the 
future  state  of  the  power  sector. 

(2)  Impact  assessment  of  C02  trading  to  business  plans’  efficiency 
and  state  licensing  policy,  within  the  context  of  Kyoto 
requirements  and  EU  directives  concerning  C02  emissions 
reduction. 

(3)  Focus  given  on  the  stability  of  the  national  grid. 

(4)  The  power  plants’  usage  factor  is  not  predetermined  but  it  may 
be  resulted  from  the  optimisation  procedure,  allowing  long 
term  planning  of  operational  intensity. 

(5)  The  Brownian  differential  of  the  stochastic  variables  are 
correlated  based  on  their  historical  data,  thus  allowing  their 
endogenous  modelling;  correlated  uncertainties  aim  to 
improve  the  performance  of  analysis  in  business  plans. 

4.  The  mathematical  model 

4.2.  Design  considerations 

The  essential  idea  of  the  model  is  the  so-called  “policy 
intervention”  in  business  and  society  [24].  Optimal  policy 
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intervention  is  attempted  through  stochastic  analysis:  the  option 
to  remain  on  existing  generation  mix  competes  with  the  option  to 
order  additional  capacities  based  on  emerging  power  production 
technologies.  An  objective  function  of  the  power  sector’s  NPV  over 
time  is  created  and  an  optimal  point  is  depicted  through  an  SQP 
optimisation  algorithm.  This  point  indicates  the  investment  entry 
timing  which  maximises  the  aggregate  system  NPV  but  also 
determines  the  upper  bound  share  of  capacities  of  each  plant 
category  (z),  allowed  to  be  ordered  (in  a  specific  time).  In  other 
words  the  state  will  be  able  to  organise  its  power  licensing 
schedule  (annually  ordered  capacities)  in  an  optimal  manner.  The 
impact  on  the  electricity  generation  industry  is  straightforward: 
let  us  suppose  that  the  variables  fluctuating  stochastically  and 
forecasted,  like  prices  of  electricity  and  fuel  (“revenues”  and 
“costs”),  are  actually  realised  in  the  future.  Also,  the  computation¬ 
ally  resolved  optimal  licensing  programmes  and/or  state  directives 
for  electricity  production  are  supposedly  violated  time-wise.  This 
means  that  some  of  the  plants  belonging  to  a  specific  plant 


theoretically  foreseen  [24].  The  entire  process  is  a  Monte-Carlo 
process  comprised  of  two  stages:  (a)  approximating  the  underlying 
multivariate  process  using  an  Euler-Marujama  numerical  solver  of 
the  vector-valued  stochastic  differential  equations  (described  in 
Appendix  B)  and  (b)  averaging  the  multiple  -  log-normally 
distributed  -  solutions;  these  outputs  feed  the  optimisation 
objective  function. 

4.3.  NPV  optimisation  and  numerical  algorithm 

As  mentioned  before,  the  model  focuses  on  the  maximisation  of 
the  system’s  NPV.  For  this  reason  an  objective  function  describing 
the  system’s  NPV  has  been  created. 

4.3 A.  Objective  function  and  variables 

The  mathematical  formulation  below  represents  the  system’s 
NPV,  considering  cumulative  annual  incomes  and  expenses 
amortised  over  operational  life-time: 


NPV(X,=1 , . . .  ,Xi=i, 0i=, 6i=l)  =  max 


Yr  I  v+leadi+li  fet 

EE  E  Pz(!  +  rr(z-v)Li.zA.z,vam,ia,,i  •  8760 

v=l  i=l  z=v+leadi 
Yr  I  v+leadi+li  fet 

-EE  E  CfuekzA  1  +  r)-{z-v)L^A,z,vam,iau  ■  8760 

v=l  i=l  z=v+leadt 
Yr  I  v+leadi+li  fet 

-EE  E  cf°mi,zO  +r)~{z~v)kz,v 

v=l  t=l  z=v+leadi 

Yr  I 

-EEcw> 

V=1  1=1 

Yr  CONVV+ leadj+ 1  i  fei  f 

-EE  E  CzAzAm+li  ■  8760  •  .(1  +  r)-{z-V  ■  pC02 

v=l  i=l  z=v+leadi  *  1 

Yr  REN  v+leadi+li  fei 


A- y  y  ^  ^i,z,v^i,z,v^m,i^i,i '  8760  •  emm[X  •  (1  +  r) 

v=l  1=1  z=v+leadi 


-( z-t ) 


P  CO, 


(1) 


category  (z)  might  operate  during  a  time  period  characterised  by 
suboptimal  conditions  due  to  the  predetermined  “revenues”  and 
“costs”  evolution  path  (i.e.  a  plant  type  (j)  of  category  (z)  which 
might  still  operate  while  its  fuel  costs  are  very  high).  Consequently, 
the  aggregate  system-wise  realised  profits  (NPV)  would  be 
suboptimal  compared  to  the  case  of  a  generic  compliance  with 
optimal  investment  timing,  strategy  and  directives,  therefore 
leaving  fewer  chances  for  profit  for  the  individual  players. 

The  maximisation  of  the  NPV  is  equivalent  to  minimisation  of 
the  costs  (C02  trading,  fuel  mix,  fixed  and  variable  costs)  while  the 
revenues  remain  unaffected,  since  they  are  derived  externally  from 
the  Monte-Carlo  simulation  of  the  stochastic  variables,  thus  giving 
an  element  of  stochastic  optimisation  to  the  model.  The  energy 
demand  is  the  critical  constraint  that  prevents  the  model  from  an 
unrealistic  infinitesimal  minimisation  of  costs,  thus  characterising 
the  model  as  a  demand-driven  study. 

4.2.  Random-walk  procedure 

The  random-walk  approximation  method  used  in  this  study 
allowed  the  simulation  of  generalised  stochastic  processes,  and 
provided  flexible  simulation  architecture.  The  procedure  sup¬ 
ported  nonlinear  relationships  which  are  commonly  found  in  SDE 
simulations.  In  the  present  study  the  Brownian  motion  differ¬ 
entials  were  correlated  based  on  historical  data  in  an  attempt  to 
simulate  the  corresponding  stochastic  variables  endogenously. 
Multiple  -  normally  distributed  -  paths  were  then  simulated  thus 
considering  the  uncertainty  of  the  model.  The  Brownian  driven 
samples  resulted  to  log-normally  distributed  output  projections  as 


The  uncertainty  is  introduced  for  the  stochastic  variables  pz, 
Cfuei,  Pco2  which  follow  the  averaged  projection  paths  generated 
through  the  Monte-Carlo  process  (described  in  Section  4.2).  The 
emmix  is  the  emissions  factor  from  the  current  conventional 
generating  mix  (in  tn  C02/MWh). 
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The  need  for  an  automated  algorithm  based  on  a  universal 
time-wise  optimisation  of  NPV  and  not  on  successive  time 
sweeps,  generated  a  triple  sum  which  feeds  the  computational 
procedure  with  the  necessary  data  that  correspond  to  the  entire 
examined  time  period.  An  innovation  of  the  model  concerns  the 
terms: 


CONV V+leadi+li  fei 


E  E  kzAzAm+u  •  8760  •  •  ■  1  •  r. 

1=1  z=v+leadi  •*  •*  * 


-(z-t) 


’  P CO, 


2534 


A.I.  Tolis  et  all  Renewable  and  Sustainable  Energy  Reviews  14  (2010)  2529-2546 


and 

RENV+leadi+li  fei 

'y  ^  ^  ^  LiZ  VOiZ  Vcimj(iii  ■  8760  •  emmix  •  (1  +  r)  ^  ^  Pco2 

i=l  z=v~\-lcudi 

which  represent  the  costs  and  revenues,  respectively,  arising  from 
the  C02  emissions  trading.  The  costs  correspond  to  the  expenses  of 
obtaining  the  required  emission  allowances  for  conventional 
power  plants,  while  the  revenues  correspond  to  the  incomes  from 
trading  the  emission  allowances  generated  by  using  renewable 
energy  sources.  The  trading  incomes  result  from  multiplication  of 
the  renewable  energy  generated  with  the  emissions  factor  of  the 
current  conventional  generating  mix  and  the  allowance  price,  i.e. 
assuming  that  the  renewable  energy  generated  replaces  energy 
that  would  otherwise  be  produced  by  the  current  conventional 
generating  mix  of  the  country. 

The  variables  of  the  objective  function  are: 


the  NPV  objective  function  and  therefore  they  are  always 
calculated  (see  Eq.  (1))  thus  reflecting  the  current  situation  in 
the  Greek  Power  Sector.  As  a  result,  the  optimisation  algorithm 
tends  to  eliminate  any  excess  energy  production,  by  iteratively 
attempting  to  reduce  the  system  costs 

4.3.2.  The  constraints  of  the  optimisation  problem 

(1)  Logical  constraints 

Usage  factor  {0^zv)  values  fluctuate  in  the  space:  0  —  @i,z,v  —  1 
( 6i  z  v  =  0  means  a  non-operational  power  plant—  0ifZtV  =  1  means  a 
maximum  possible  operational  time).  Also  the  annual  orders  for 
new  power  plants  Xiv  should  be  positive  numbers:  Xiv  >  0,  V i,  v. 

(2)  Natural  resource  availability 


(i)  the  orders  (capacities)  realised  in  the  year  ( v ):  Xiv  (MW)  for 
investing  on  new  power  plants  of  a  specific  technology  and/or 
fuel  type  (i)  and, 

(ii)  the  usage  factor  0iz  y  of  the  different  types  of  installed  power 
plants  (i)  operating  during  year  (z)  and  ordered  during  year  (v). 


It  is  noted  that  0iz  v  is  requested  only  for  the  power  generating 
units  fired  by  expensive  fuels  (e.g.  lignite,  natural  gas),  otherwise 
(e.g.  solar  PV,  wind-farms)  it  is  imposed  0itZ<v  =  1.  This  means  that 
the  power  plants  with  zero  fuel  cost  should  be  used  whenever  their 
energy  source  is  available,  thus  reflecting  real  world  practice  while 
slightly  reducing  the  variables  and  the  computational  cost. 

The  investment  costs  CIi  v,  of  a  power  plant  type  (i)  in  the  year 
(v)  depend  on  technical  advances  arising  from  long  periods  of 
cumulative  experience  on  construction  of  such  power  production 
units.  This  fact  is  expressed  by  the  learning  curves  methodology 
which  can  be  mathematically  formulated  as  follows: 


CIi,v  —  C/j52000 
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The  investment  expenses  of  year  2000  are  used  as  a  reference 
value. 

Two  scenarios  were  considered,  concerning  the  electricity 
demand: 


-  In  the  first  scenario  (forecasted  prices  and  demand  or  briefly  FPD) 
the  future  loads  are  forecasted  using  historical  data  (the  relative 
numerical  approach  is  described  in  Section  4.2). 

-  In  the  second  (elasticity  related  prices  and  demand  or  briefly 
ELPD)  the  demand  is  related  to  its  price  with  the  following 
elasticity  formula  (similar  to  Ref.  [2]): 


dz  =  a ■  pzE  (3) 

This  formula  expresses  a  condition  realised  by  the  assumption: 
“Electricity  demand  depends  only  on  its  price  and  does  not  have 
any  antagonists  or  counterparts  impacting  its  price.  Also,  its 
elasticity  (E)  may  be  considered  as  constant”.  Therefore  in  both 
scenarios,  the  total  demand  and  the  forward  prices  are  related  in  an 
attempt  to  capture  market  dynamics  endogenously. 

Excessive  produced  energy  considerations :  It  is  noted  that  energy 
produced  in  excess  to  the  demand  is  a  condition  occasionally 
occurred  during  the  numerical  iterations  of  the  optimisation 
algorithm  until  it  converges  to  an  optimal  solution.  Any  excess 
energy  produced  is  assumed  to  have  zero  value,  as  it  cannot  be 
exploited  by  the  system.  On  the  other  hand,  the  production  costs  of 
the  non-served  energy,  contribute  (being  a  part  of  the  expenses)  in 


Some  technologies  of  power  production  are  constrained  by 
limited  natural  resources  and/or  fuel  reserves  which  determine  the 
upper  limit  of  the  electricity  generation  potential.  The  wind 
generators,  the  hydro  plants  and  the  geothermal  plants  are 
indicative  examples  of  such  types  of  electricity  generation  whose 
resources  present  an  upper  potential  limit.  This  constraint  can  be 
represented  mathematically  as  follows: 

z+leadj+li  fej 

Y  XitV  <  dyrij  V  i  (4) 

V=z+leadi 

(3)  Kyoto  protocol  and  EU  Directive  2001/77/EC  constraints 


Since  the  emissions  trading  has  been  engaged  in  the  model’s 
objective  function,  additional  constraints  should  be  entered 
concerning  the  compliance  with  the  Kyoto  protocol  requirements 
and  the  EU  Directive  2001 /77/EC.  These  constraints  are  obligatory 
for  every  participating  country.  Greece  is  required  to  meet  the 
following  targets: 


(3a)  the  share  of  renewable  energy  sources  should  be  at  least  20,1  % 
of  the  total  electricity  production  by  the  year  2010.  In 
mathematical  formulation: 
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(3b)  the  share  of  renewable  energy  sources  should  be  at  least  30% 
of  the  total  electricity  production  by  the  year  2020.  In 
mathematical  formulation: 
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(4)  Social  constraints 

Meeting  the  demand  target  is  a  basic  condition  ensuring  social 
acceptance  of  power  production  at  a  national  level.  In  the  present 
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study  the  total  annual  electricity  production  is  imposed  to  be  at 
least  as  much  as  the  predicted  annual  electricity  demand.  This 
constraint  is  split  in  two  discrete  branches: 

(4a)  Total  demand  target  (energy):  As  stated  before,  two  discrete 
scenarios  (FPD  and  ELPD)  have  been  considered  in  the  model 
regarding  the  projection  of  demand.  In  both  scenarios,  the 
electricity  generation  (in  terms  of  energy)  should  exceed  the 
projected  demand  thus  being  able  to  serve  the  yearly  demand 
aggregations: 
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The  rationale  for  this  constraint  was  based  on  the  assumption 
that  a  possible  generic  failure  of  power  supply  (system  black-out) 
would  result  to  an  excessive  social  cost.  Therefore  this  case  should 
be  avoided  in  the  expense  of  additional  capacity  orders  causing  a 
system  NPV  reduction.  An  additional  safety  factor  has  been 
implemented,  namely  the  forecast  error.  In  the  FPD  scenario,  the 
MAD  (mean  absolute  deviation)  indicator  has  been  used  as  a 
measure  of  success  of  the  forecast  for  a  predetermined  time  period 
used  for  its  validation.  However,  -  given  the  increased  uncertainly 
of  the  lognormal  forecasted  demand  arrays  -,  this  MAD  value  has 
been  scaled  up  in  order  to  adjust  to  the  entire  time  horizon  of  the 
model.  In  the  ELPD  scenario  a  value  derived  from  the  elasticity 
formula  and  the  electricity  price  forecast  MAD  has  been  used. 

(4b)  Peak-power  demand  target:  The  constraint  implies  that  the 
aggregate  installed  capacities  should  be  able  to  serve  the 
projected  peak-power  demand.  The  peak-power  demand  - 
annual  past  observations  -  follows  a  linear  low-slope  trend 
and  therefore  it  has  been  linearly  projected  to  the  future.  The 
availability  and  usage  factors  (except  those  of  the  wind 
turbines  and  the  solar  PVs)  have  been  considered  to  be  equal 
to  1  meaning  full  load  operating  conditions  with  available 
energy  potential  for  the  peak-power  demand  period.  A  safety 
factor  (generation  reserve  margin  taking  values  in  the  range 
0,1 -0,2)  has  been  further  multiplied  to  the  peak-power 
demand  in  order  to  account  for  a  possible  non-availability  of 
wind  or  solar  energy  or  the  possibility  of  some  plants’  failure: 
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Two  discrete  assumptions  have  been  formed  for  each  one  of  the 
demand  scenarios  concerning  the  peak-power  demand  estima¬ 
tions: 

(i)  FPD  scenario:  The  peak-power  demand  -  as  observed  in  each 
past  year  -  followed  a  linear  trend  and  therefore  has  been 
linearly  projected  to  the  future. 

(ii)  ELPD  scenario:  This  scenario  generally  predisposed  a  moderate 
growth  concerning  the  demand  uncertainty.  Therefore,  the 
peak-power  demand  has  been  related  to  the  electricity  prices 
in  the  same  manner  with  the  total  electricity  demand  forecast 
(using  Eq.  (3)— elasticity  related  prices  and  demand). 

This  constraint  might  be  able  to  model  any  possible  case  of  the 
system’s  reliability  limits  (either  total  energy  or  peak  power 
depending  on  which  constraint  proves  to  be  stricter).  The 
generation  reserve  margin  concept  provides  additional  safety  on 
the  model  and  has  been  also  adopted  in  the  works  [2,5]. 

(5)  Stability  of  the  national  grid 


It  is  well  known  that  some  power  generation  technologies 
which  depend  strongly  on  natural  resources  cannot  be  the  base¬ 
load  workhorses  of  electricity  production.  Such  technologies  are: 
the  wind  turbines,  the  solar  PVs,  the  geothermic  plants  and  in 
general  most  of  the  types  of  Power  Generators  that  rely  on 
renewable  resources.  Despite  of  their  short  setup  periods  and 
cheap-or-no  need  for  fuels,  they  often  suffer  from  resource 
unavailability  and  unpredictability  and  can  severely  affect  the 
national  grid’s  stability.  In  the  present  study,  it  has  been  assumed 
that  the  cumulative  annual  electricity  energy  production  based  on 
such  technologies  should  be  less  than  50%  of  the  total  annual 
demand: 
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(6)  The  limits  of  state  licensing  procedure 


The  upper  bound  of  the  annually  ordered  capacities  from  each 
one  technology  should  be  less  than  600  MWei,  or  Xiv  <  600  Vi,  v. 
This  value  was  a  conclusion  mainly  based  on  the  past  experience 
for  Greece,  reflecting  the  limits  of  its  economy’s  robustness  in 
conjunction  with  power  demand. 

5.  The  computational  procedure 

The  data  structures  have  been  modified  in  order  to  create  arrays 
with  compatible  time  discretisation.  Initially,  the  past  data  were 
available  in  different  time  intervals.  Electricity  prices  (Euro / 
MWhei)  and  demand  (MW)  were  available  in  an  hourly  basis;  C02 
prices  were  available  in  a  daily  basis,  while  fuel  prices  and  central 
bank  rates  were  available  in  a  monthly  basis.  This  was 
contradictory  with  the  requirement  of  supplying  the  optimisation 
algorithm  with  annual  inputs.  Compatible  data  structures  were 
mandatory  in  order  for  the  optimisation  code  to  be  usable. 
Therefore,  prior  to  the  optimisation  algorithm  runs,  the  historical 
data  were  subject  to  the  following  modifications: 

(a)  Initially,  some  fuel  prices,  electricity  prices  and  demand  were 
missing.  Linear  interpolation  techniques  have  been  applied  to 
overcome  this  problem. 

(b)  The  granularity  selected  for  creating  the  history  of  stochastic 
variables  was  based  on  the  widest  available  data  sets  [35] 
which  for  the  available  data  were  one  month-periods.  This 
means  that  e.g.  daily  prices  can  be  averaged  in  monthly  prices, 
but  monthly  prices  could  not  be  transformed  in  daily  prices 
unless  a  regression  or  interpolation  technique  is  applied  in  the 
expense  of  inaccurate  variations  modelling.  Therefore  the 
historical  data  were  averaged  in  monthly  basis  unless  they 
were  already  provided  in  that  granularity.  Concerning  a 
forward  electricity  contract  with  maturity  T1  and  delivery 
period  [Tl,  T2],  an  approximation  of  monthly  averages  F(T) 
(forward  prices)  can  be  realised  if  daily  (hourly)  prices  /(t)  are 
initially  available  within  daily  (or  hourly)  intervals  dt  according 
to  the  formula: 

F(T)  =  T2[n  fd)dt  Pe[Tl,T2]  (10) 

It  is  noted  that  while  the  electricity  prices  were  actually 
averaged,  this  was  not  the  case  for  the  demand  loads.  Demand 
forecasts  were  derived  by  multiplying  the  demand  observa¬ 
tions  d(t)  by  their  occurrence  time  interval  dt,  and  aggregating 
the  products  within  each  month  thus  resulting  to  total  monthly 
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Table  1 

List  of  technology  related  parameters  (source:  Greek  Ministry  of  Development). 


Integrated 

coal 

gasification 

Oil 

Natural 

gas 

Lignite 

Biomass 
gasifier 
STAG  (NH) 

Solar 

PV 

Large 
onshore 
wind  turbines 

Medium 
and  high 
head  hydro 

Hydro 

pumped 

storage 

Geothermal 

Investment  cost  (103€/MWei)  of  the  year  2000 

1376 

1127 

1062 

1888 

2674 

6555 

1245 

2491 

3736 

1350 

Fixed  cost  of  year  2000  (operational, 
maintenance,  insurance,  etc.)  (in  103€/MWei) 

70,80 

31,47 

27,53 

48,51 

262,2 

26,88 

23,60 

24,91 

49,82 

34,09 

Availability  factor 

0,75 

0,75 

0,75 

0,82 

0,75 

0,90 

0,90 

0,85 

0,92 

0,70 

Capacity  factor 

0,80 

0,80 

0,65 

0,75 

0,80 

0,15 

0,25 

0,34 

0,40 

0,90 

Learning  rate 

0,05 

0,01 

0,01 

0,01 

0,15 

0,20 

0,10 

0,00 

0,00 

0,00 

Commissioning  time 

4 

3 

3 

4 

3 

2 

1 

10 

10 

3 

Efficiency  factor 

0,45 

0,38 

0,56 

0,36 

0,29 

1 

1 

1 

1 

1 

Fuel  C02  emissions  (tn  C02/MWh  fuel) 

0,354 

0,278 

0,202 

0,363 

0 

0 

0 

0 

0 

0 

Operational  life-time  (years) 

30 

30 

30 

40 

20 

30 

30 

45 

45 

20 

energy  demand  estimations  (D(T)  in  MWhei): 

rT2 

D(T)  =  /  d(t)dt  Te  [71,72]  (11) 

(c)  When  in  monthly  form,  the  variables  were  mathematically 
correlated  using  correlation  factors  derived  by  their  past 
monthly  history  (either  averaged  or  original)  in  a  similar  vein 
to  [36].  Actually,  the  Brownian  differentials  of  fuel,  electricity 
prices  and  electricity  demand  were  correlated  (mutually  as 
well  as  with  the  corresponding  Brownian  differentials  of  the 
interest  and  inflation  rates).  Especially  concerning  the  histori¬ 
cal  C02  allowance  prices,  they  have  been  correlated  with  the 
historical  electricity  prices  as  has  been  adopted  in  Refs.  [7,8]. 
The  statistical  parameters  required  for  the  GBM  (or  MR) 
simulation  (drift,  standard  deviation,  correlation  factors)  were 
extracted  by  the  monthly  arrays  of  the  historical  data. 

(d)  The  GBM  (or  MR)  simulation  of  the  stochastic  variables 
produced  multiple  forecasts  through  the  Monte-Carlo  algo¬ 
rithm  described  in  Section  4.2.  These  projections  were  then 
averaged  in  a  yearly  basis.  Monthly  or  even  shorter  granula¬ 
rities  have  been  also  adopted  in  the  past  (e.g.  [21]).  Yearly 
discretised  studies  have  been  also  conducted  by  [2,5]. 

The  computational  procedure  used  for  the  problem’s  numerical 

modelling  can  be  summarised  to  the  following  steps: 

Step  1.  Determination  of  problem’s  inputs  and  technological 
details,  such  as:  setup  periods,  lead  time  (construction), 
operational  life-time,  basic  investment  and  operational 
costs,  and  load  factors. 

Step  2.  Simultaneous  prediction  of  electricity  selling  price,  elec¬ 
tricity  demand,  C02  allowance  and  fuel  prices.  Correlated 
Brownian  differentials  -  based  on  appropriately  structured 
past  data  -  are  used  through  a  multivariate  Monte-Carlo 
process  in  which  an  Euler-Marujama  solver  is  used  to 
produce  multiple  solutions  of  either  GBM  or  MR  modelling. 

Step  3.  Supplying  the  optimisation  algorithm  with  the  numerical 
forecasts  generated  in  step  2  (after  appropriate  post¬ 
processing  of  data  structures). 

Step  4.  Optimisation  of  the  future  power  sector  structure  using  an 
iterative  SQP  method.  Reconstitution  of  future  electricity 
production  by  plant  type  (and  consequently  by  fuel  source). 


5.2.  The  model’s  input  data 

Three  basic  parameters  characterise  the  entire  model:  time, 
inflation  and  interest  rates.  The  time  period  of  the  study  is  50  years 
but  only  the  first  25  years  might  include  valid  results  as  will  be 
explained  later  on.  The  inflation  rate  is  assumed  to  be  4%,  close  to 


the  current  mean  average  of  Greece,  while  the  interest  rate,  which 
includes  the  inflation  and  a  technological  risk,  is  assumed  to  be  5%. 
Another  set  of  inputs  necessary  for  the  model  was  the  initial  status 
of  the  Greek  Electricity  Generation  system  as  of  the  year  2001  [10]. 
The  entire  field  of  technology  related  inputs  needed  for  the 
iterative  computational  procedure  is  presented  in  table  form 
(Table  1),  including  information  for  the  various  parameters  and 
variables. 

5.2.  Prediction  of  electricity  prices  and  demand 

Among  the  various  possible  combinations  of  input  parameters 
and  assumptions  tested  with  the  model’s  code,  only  two  scenarios 
were  selected  for  presentation.  These  two  case  studies  represent 
reasonable  hypotheses  for  the  future  of  deregulated  electricity 
markets: 

5.2.2.  Case  study  A— explicitly  forecasted  electricity  prices  and 
demand  (FPD) 

In  deregulated  markets  the  competition  for  domination  usually 
pushes  goods’  prices  down.  Electricity  -  being  strongly  inelastic  - 
could  be  influenced  by  prices’  fluctuation  only  in  rare  circum¬ 
stances.  Its  evolution  pattern  usually  follows  the  actual  social 
needs  through  time.  For  this  reason,  the  first  case  study  was  based 
on  the  assumption  of  increasing  demand  drift  (slope),  as  observed 
in  the  historical  data  (step  1)  of  demand  acquired  by  the  Hellenic 
Transmission  System’s  Operator  [37].  The  electricity  prices 
(national  grid’s  marginal  prices)  were  also  based  on  historical 
data  acquired  by  the  same  source.  This  data  actually  constitutes 
the  payments  realised  by  HTSO  S.A.  to  the  electricity  producers. 
The  historical  data  used  for  both  variables  was  initially  available  in 
an  hourly  basis  for  years  2001-2008,  but  they  were  computation¬ 
ally  modified  as  described  in  the  introduction  of  chapter  5.  A 
Monte-Carlo  procedure  (multiple  solutions  of  a  GBM  model  as 
described  in  Section  4.2)  produced  the  annual  estimation  of  total 
demand  loads  and  forward  electricity  prices  which  can  be  seen  in 
Figs.  1  and  2,  respectively. 

The  standard  deviation  for  the  historical  data  of  peak  electricity 
demand  has  been  calculated  equal  to  0,037  while  its  mean  drift 
was  0,018. 

The  stochastic  differential  standard  deviation  (defined  by 
dividing  standard  deviation  with  mean  stochastic  differential) 
for  the  historical  electricity  prices  has  been  calculated  equal  to 
14,97. 

5.2.2.  Case  study  B— elasticity  related  horizon  of  electricity  prices  and 
demand  (ELPD) 

The  national  grid  marginal  prices  were  extracted  by  HTSO  S.A. 
as  in  the  first  case  study.  The  same  GBM  model  and  the  same 
Monte-Carlo  numerical  procedure  were  used  in  order  to  produce 
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Fig.  1.  The  forecasted  electricity  demand  for  FPD. 
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Fig.  3.  The  forecasted  electricity  demand  for  ELPD. 


their  future  field.  The  nominal  electricity  prices  obtained  were  then 
transformed  to  real  prices  (constant  prices  with  reference  to  the 
year  2001),  based  on  the  argument  which  implies  that  during  an 
inflation  period  the  expenses  rise  with  almost  the  same  rate  as  the 
incomes  do,  so  prices  of  goods  should  be  referenced  to  a  constant 
time,  in  order  to  reveal  their  true  future  influence  to  demand, 
through  elasticity.  The  prices’  reference  to  a  constant  time  could  be 
solely  carried  out  using  the  assumed  inflation  rate,  unless  an 
investment  cash-flow  had  been  considered  which  would  require 
interest  rate-based  transformation. 

A  constant  elasticity  equal  to  0,4  was  assumed  for  the  entire 
studied  time  period.  The  elasticity  equation  (3)  was  then  applied  in 
order  to  project  the  future  electricity  demand.  It  is  noted  that  in 
order  to  compute  the  elasticity  formula’s  base,  a  simple  algebraic 
equation  was  solved  that  related  the  price  and  the  demand  of 
electricity  (both  known  values)  during  the  first  year  of  examina¬ 
tion.  The  assumptions  made  above  were  based  on  the  initial 
hypothesis  of  non-existence  of  electricity  substitutes  impacting  its 
price.  This  method  finally  resulted  to  a  moderate  increase  of  total 
demand  as  can  be  seen  in  Fig.  3.  The  forecasted  forward  prices 
retain  the  same  path  as  in  the  previous  case  (FPD).  Average 
electricity  prices  have  been  used  for  deriving  the  total  demand 
(including  peak  loads). 

In  both  cases  the  expected  demand  increases  with  time  but  its 
slope  is  steeper  in  the  FPD  case.  In  the  ELPD  case  the  demanded 
loads  are  moderately  increasing  as  a  result  of  the  constant 
elasticity  assumption.  Both  cases  comply  with  a  reasonable 
hypothesis:  “electricity  demand  keeps  up  with  economies’ 
robustness  and  financial  growth”  [4,38].  It  is  noted  that  Greece 


experienced  for  almost  a  decade  a  constant  economic  growth  in  the 
levels  of  3%  or  higher  [39]. 

5.3.  C02  allowance  prices 

The  C02  allowance  prices  contribute  to  the  system’s  expenses 
and  revenues  as  seen  in  Eq.  (1).  These  prices  present  a  constant  - 
but  slow  -  reverting  to  a  long  run  mean  derived  by  historical  data 
[40],  As  a  consequence,  the  forecasting  of  these  prices  was 
performed  by  simulating  an  (MR)  type  of  SDE  (Eq.  (3*)  of 
Appendix  B).  The  statistical  analysis  indicated  a  medium  correla¬ 
tion  with  electricity  prices  (correlation  factor  0,62).  The  forecasted 
C02  allowance  prices  can  be  seen  in  Fig.  4. 

5.4.  Fuel  prices 

Fuel  prices  contribute  significantly  to  the  finances  of  a  power 
plant.  Therefore,  they  should  be  forecasted  paying  particular 
attention  on  the  correlation  of  their  evolution  patterns.  This  was 
done  by  simulating  the  GBM  type  SDEs  (Eq.  (2*)  of  Appendix  B).  The 
solution  was  repeated  for  almost  500  trials  within  the  frame  of  an 
iterative  Monte-Carlo  simulation  using  correlated  stochastic 
differentials  as  explained  in  the  description  of  the  mathematical 
model. 

The  historical  data  of  fuel  prices  cover  the  time  period  2001- 
2008  and  provide  the  necessary  information  about  the  statistical 
parameters  required  by  the  GBM  simulation  (mean  drift,  standard 
deviation,  correlation  factors,  etc.).  Past  data  sources  were  among 
others:  the  Hellenic  Ministry  of  Development  (hard  copies)  as  well 
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Fig.  2.  The  forecasted  forward  electricity  prices. 
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Fig.  4.  The  forecasted  C02  allowance  prices. 


Table  2 

Stochastic  differentials’  correlation  coefficient  of  fuel  prices  extracted  from  historical  values. 


Hardcoal 

Oil 

Natural  gas 

Lignite 

Biomass 

Geothermal 

Hydroelectric 

Hardcoal 

1 

0,87 

0,89 

-0,11 

-0,57 

1,00 

0,78 

Oil 

0,87 

1 

0,67 

-0,20 

-0,30 

0,88 

0,90 

Natural  Gas 

0,89 

0,67 

1 

-0,39 

-0,88 

0,89 

0,48 

Lignite 

-0,11 

-0,20 

-0,39 

1 

0,57 

-0,13 

0,13 

Biomass 

-0,57 

-0,30 

-0,88 

0,57 

1 

-0,57 

-0,04 

Geothermal 

1,00 

0,88 

0,89 

-0,13 

-0,57 

1 

0,78 

Hydroelectric 

0,78 

0,90 

0,48 

0,13 

-0,04 

0,78 

1 

as  Ref.  [41].  Contractual  lignite  prices  were  available  and  acquired 
from  private  mines.  Concerning  the  biomass  price,  the  technique 
established  in  Ref.  [42]  was  applied,  which  is  a  Reverse  Logistics 
approach  augmented  by  a  Holistic  Activity  Based  Costing  proce¬ 
dure.  Having  as  a  starting  point  the  price  calculated  with  this 
technique  for  the  year  2001 ,  a  mean  drift  equal  to  the  inflation  rate 
was  further  applied  in  order  to  produce  the  biomass  price 
evolution  till  the  year  2008  (last  year  of  input  data).  The  table 
of  stochastic  differential  correlations  derived  by  the  previously 
analysed  historical  data  of  fuel  prices  is  presented  (Table  2). 

The  resulting  fuel  price  evolution  is  presented  in  the  following 
graph  (Fig.  5)  as  an  outcome  of  the  GBM  SDE  solution  augmented 
by  the  Monte-Carlo  simulation.  It  is  emphasised  that  the  fuel  prices 
include  their  mining  and/or  harvesting  costs,  their  logistical  costs 
(transportation,  handling  and  warehousing)  as  well  as  the  costs  for 
their  production  or  processing  to  their  final  form. 

The  prices  presented  in  the  above  diagram  (values  given  in  €/ 
MWhei)  follow  in  general  their  historically  observed  trend.  For 


example  oil  prices  are  slowly  reverting  to  a  constant  equilibrium 
which  is  a  typical  behaviour  for  oil  products  [24].  The  anticipated 
coal  prices  are  decreasing,  which  is  in  line  with  studies  concerning 
the  global  coal  market  [43].  On  the  other  hand,  anticipated  lignite 
prices  present  an  increasing  trend.  The  gradually  depleting  Greek 
reserves,  might  actually  lead  to  harder  winning  efforts,  thus 
justifying  this  behaviour.  The  statistical  analysis  indicated  a 
definite  correlation  of  lignite  and  electricity  prices  (0,77).  Gas 
and  electricity  prices  presented  a  slightly  weaker  correlation  factor 
0,71.  The  previously  analysed  and  processed  data  are  gathered  as 
inputs  for  the  optimisation  algorithm. 

5.5.  The  time  periods  of  the  study 

Within  the  frame  of  this  study  one  could  easily  establish  three 
discrete  time  periods: 

The  initial  aim  of  the  study  was  to  calculate  the  optimum 
capacity  orders  for  the  next  25  years,  meaning  till  the  year  2033. 
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Fig.  5.  The  forecasted  fuel  prices. 
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producing  technology,  in  order  to  extract  more  accurate  results 
one  should  modify  the  inputs  by  selecting  a  much  wider  time 
period  (i.e.  till  the  year  2080),  thus  being  able  to  investigate  the 
influence  of  power  plants  to  be  built  in  the  distant  future.  However, 
adopting  an  even  wider  timeframe  would  lead  to  two  main 
obstacles:  (i)  In  its  current  form  the  computational  code  examines 
the  50-year  period  and  its  CPU  runtime  approaches  the  limits  of 
the  acceptable,  (ii)  A  wider  time  frame  would  actually  cancel  the 
validity  of  the  random-walk  forecasting  procedure.  Since  its  nature 
implies  a  constantly  increasing  uncertainty  of  the  log-normally 
distributed  results  [24],  it  can  be  easily  understood  why  the  time 
limits  of  implementation  should  be  tight  (Fig.  6).  Therefore, 
although  the  code  runs  for  50  years,  only  the  results  concerning  the 
time  period  2009-2033  may  be  considered  valid,  as  was  initially 
planned. 

6.  Results  of  the  code 


However,  the  iterative  optimisation  procedure  should  run  for  a 
wider  timeframe  because  some  technologies  require  significant 
construction  times  (up  to  10  years).  Therefore,  power  loads  and 
their  financial  influence  should  be  calculated  at  least  till  the  year 
2043.  Another  set  of  7  years  were  included  in  order  to  let  the 
computational  code  eliminate  the  time-boundary  effects.  For  this 
reason,  the  upper  time  limit  was  imposed  to  reach  the  year  2050. 
Despite  of  true  operational  life-times,  used  for  each  power 


The  results  of  the  model  are  presented  in  the  following  sub¬ 
sections. 

6.2.  The  FPD  case 

The  capacity  additions  appear  to  have  a  complex  profile  as  a 
function  of  time.  In  the  valid  range  of  results  there  is  no  plant 
category  expected  to  dominate  compared  to  others.  Almost  all 

■  Geothermal  power  plant 
Hydro  pumped  storage 
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Large  onshore  wind  turbine 
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Fig.  7.  The  calculated  capacity  orders  for  FPD. 
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Fig.  8.  The  calculated  usage  factor  for  FPD. 
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Fig.  9.  The  cumulative  energy  produced  (FPD)  in  comparison  with  forecasted  demand. 


Fig.  10.  The  expected  C02  emissions  per  produced  MWhei  for  FPD. 


plant  types  should  be  ordered  to  contribute  in  the  increased  power 
demands  assumed  for  the  FPD  case.  The  expected  yearly  capacity 
orders  can  be  seen  in  Fig.  7. 

The  expected  usage  factor  smoothly  oscillates  in  high  levels  for 
the  same  reason  (Fig.  8). 

In  Fig.  9  the  cumulative  yearly  electricity  production  is  shown. 
It  is  interesting  to  note  that  electricity  production  succeeds  to  meet 
the  demand  target  shown  in  the  upper  edge  of  the  curves  with  a 
line. 


An  interesting  decrease  of  C02  emissions  per  MWhei  produced 
is  expected  and  this  can  be  seen  in  the  graph  of  Fig.  10.  The 
anticipated  increasing  penetration  of  renewable  energy  sources  - 
shown  in  Fig.  1 1  -  might  justify  this  behaviour.  This  is  also  an 
indication  that  the  optimisation  procedure  attempted  to  lower  the 
emission  costs  thus  improving  the  system’s  economy  (i.e. 
compared  to  the  current  situation  of  year  2009). 

The  computational  code  succeeded  to  follow  the  required 
increase  of  renewable  energy  penetration  against  the  conventional 
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Conventional 

Renewables 


Fig.  11.  The  expected  penetration  of  renewable  energy  (including  hydroelectricity)  in  the  Greek  Power  Sector. 
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power  plants.  More  specifically,  starting  from  the  levels  of  2001 
(close  to  8%)  the  renewables’  contribution  rises  to  the  current 
levels  (close  to  1 5%  in  2008),  and  also  manages  to  comply  with  the 
EU  directives  (20,1%  in  2010  and  30%  in  2020). 

The  calculated  costs  breakdown  is  presented  in  Fig.  12.  The 
expected  share  of  C02  expenses  -  due  to  emissions  trading  -  is 
moderate  and  slightly  rising  with  time.  However  its  percentage  in 
comparison  with  the  expected  fuel  costs  (almost  half)  is  significant 
enough  to  impact  the  oncoming  investment  decisions. 

The  expected  financial  balance  is  presented  in  Fig.  13  for  the 
entire  system.  The  NPV  of  the  balance  curve  corresponds  to 
approximately  5  x  109  €  as  can  be  seen  in  Fig.  20. 


6.2.  The  ELPD  case 

The  ELPD  case  produced  a  different  decomposition  of  the  Greek 
electricity  sector.  Within  the  range  of  valid  forecasts  (years  2009- 
2033,  Fig.  14)  hydroelectric  and  large  onshore  wind  turbines 
should  dominate  the  new  orders,  to  maintain  the  medium  levels  of 
electricity  demand  arising  from  the  constant  elasticity  hypothesis. 
Solar  PV  plants  should  be  installed  in  the  near  future  whereas  the 
main  base-load  technology  will  remain  the  existing  lignite  plants. 

As  can  be  seen  in  Fig.  15  the  usage  factor  should  vary  in  lower 
levels  comparing  to  those  observed  in  the  FPD  case  especially  for 
the  lignite  fired  plants.  This  happens  mainly  due  to  the  increasing 


Fig.  13.  The  expected  system’s  financial  balancing  for  FPD. 
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Fig.  14.  The  calculated  capacity  orders  for  ELPD. 


costs  forecasted  for  the  case  of  lignite  and  the  large  share  of  lignite 
power  plants  at  the  starting  period. 

High  90%  are  optimal  levels  for  usage  factors  for  all  but  lignite 
fired  plants  which  should  operate  at  80%  or  lower  for  most  of  the 
time.  The  slight  electricity  demand  increase  of  ELPD  case  is  shown 
in  Fig.  16  in  comparison  with  the  cumulative  electricity  produced: 
an  obviously  easier  target  to  meet  compared  to  the  steep  electricity 
demand  increase  of  the  FPD  case. 

Due  to  the  fact  that  the  required  additions  are  minor,  the 
penetration  of  renewable  energy  sources  is  slower  than  in  the  FPD 


case.  However  the  expected  projection  of  emissions  remains  in 
significantly  low  levels  too.  As  a  result  the  expected  C02  emissions / 
MWhei  produced  are  decreasing,  as  shown  in  Fig.  17  thus 
indicating  a  consequent  economy  improvement  due  to  the 
optimisation  algorithm. 

The  expected  expenses  and  financial  balancing  breakdowns  are 
shown  in  Figs.  18  and  19,  respectively.  The  C02  trading  costs  still 
contribute  significantly  to  the  system  costs. 

It  has  to  be  noted  that  in  Figs.  18  and  19,  current  (nominal) 
prices  are  presented  meaning  that  their  present  value  is 
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Fig.  16.  The  cumulative  energy  produced  (ELPD)  in  comparison  with  forecasted  demand. 
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Fig.  17.  The  expected  C02  emissions  per  produced  MWhei  for  ELPD. 


significantly  lower,  thus  resulting  to  a  final  NPV  close  to  7  x  109  € 
as  shown  in  Fig.  20. 

6.3.  Sensitivity  analysis 

Despite  of  the  model’s  robustness  for  most  of  the  input 
parameters,  the  system’s  aggregate  NPV  could  be  severely 
influenced  by  the  assumption  of  the  interest  rates.  Therefore, 
the  model’s  sensitivity  should  be  further  investigated  for  a  wider 
range  of  rates  assumed.  The  computational  procedure  has  been 


successively  performed  for  acceptable  rates  in  the  range  4-8%  with 
increasing  step  equal  to  1%.  The  resulting  sensitivity  is  shown  in 
Fig.  20  for  both  scenarios. 

Lower  interest  rate  values  lead  to  higher  NPVs  as  reasonably 
expected  while  a  more  moderate  decrease  appears  for  interest 
rates  exceeding  6%.  However,  negative  aggregate  NPV  is  noted  for 
high  interest  rate  values,  alerting  for  deficit  systemic  operation.  It 
is  noted  also  that  the  ELPD  scenario  proves  to  be  less  sensitive 
compared  to  FPD  mainly  due  to  the  significantly  lower  need  for 
future  investments  and  fuel  cost. 


Fig.  18.  The  projected  annual  costs  and  their  expected  breakdown  for  ELPD. 
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Fig.  19.  The  expected  system’s  financial  balancing  for  ELPD. 


interest  rate  (%) 

Fig.  20.  Sensitivity  of  the  model  for  the  ELPD  and  FPD  cases. 


7.  Conclusions 

A  computational  model  has  been  created  in  order  to  calculate 
the  optimal  power  generation  capacity  additions  in  the  context  of 
future  uncertainties.  The  model  incorporates  the  calculation  of  the 
optimal  usage  factor  of  the  power  plant  as  well  as  the  impact  of 
emissions  trading  to  the  finance  of  the  system.  The  aggregate 
system  NPV  constitutes  the  performance  criterion  of  the  study. 
Focus  has  been  shown  on  the  electricity  sector  of  Greece.  Two 
discrete  scenarios  have  been  examined  in  detail.  In  the  first  one, 
the  electricity  demand  was  forecasted  directly  from  historical  data 
(FPD)  while  in  the  second,  a  predetermined  elasticity  defined  the 
future  path  of  the  electricity  demand  (ELPD). 

In  the  FPD  scenario,  the  anticipated  generation  mix  consists  of 
existing  but  proven  technologies  as  well  as  new  emerging 
technologies  based  on  renewable  energy  sources.  The  steep 
increasing  of  the  electricity  demand  imposes  the  requirement 
for  power  generation  originating  from  all  the  available  technolo¬ 
gies,  but  not  uniformly  distributed.  The  low-or-zero  fixed  and  fuel 
cost  of  wind  generators  and  hydroelectric  plants  leads  to  a 
significant  percentage  of  the  required  capacities  being  based  on 


these  technologies.  However,  fossil-fueled  power  plants  like  lignite 
and  natural  gas  fired  plants,  should  retain  their  position  as 
strategic  base-load  fuels,  as  is  the  current  situation  in  Greece.  The 
remaining  renewable  source  technologies  like  solar  PVs  and 
biomass-fired  plants  are  challenged  to  this  race  and  they  may  not 
gain  a  higher  part  of  power  generation  unless  experience  acquired 
in  the  future  reduces  their  investment  cost.  The  usage  of  most  plant 
types  should  remain  in  high  levels  close  to  90%  in  order  to  meet  the 
demand  target.  The  produced  electricity  -  as  calculated  by  the 
model  -,  could  comply  with  the  Kyoto  protocol  requirements, 
concerning  their  future  penetration  against  conventional  power 
plants:  a  share  exceeding  20,1%  and  30%  of  the  total  production  for 
the  years  2010  and  2020,  respectively  should  be  based  on 
renewable  energy  sources.  Additionally,  the  stability  of  the  system 
is  ensured  due  to  the  restriction  of  renewable  energies  far  below 
the  limit  of  50%  of  electricity  demand. 

In  the  ELPD  scenario  the  demand  target  is  moderate  and 
therefore  easier  to  be  met  by  the  optimisation  algorithm.  The  low 
investment  and  fuel  cost  of  hydroelectric  plants  favours  their 
future  exploitation.  Additional  plant  types  based  on  renewable 
sources  of  energy  like  Solar  PVs  and  biomass-fired  plants  would  be 
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needed  in  low  percentages,  since  the  demand  increase  is  not  steep. 
For  the  same  reason,  the  usage  factor  of  future  and  existing  lignite 
plant  types  should  be  moderate  and  rather  lower  than  those 
calculated  in  the  FPD  scenario. 

The  model  proved  to  be  robust  in  most  of  the  input 
parameters  but  sensitive  on  the  imposed  interest  rates.  Regard¬ 
ing  the  performance  of  the  two  investigated  scenarios,  the  ELPD 
scenario  appears  to  be  more  efficient  in  terms  of  the  aggregate 
anticipated  NPV.  The  moderate  requirements  for  electricity 
demand  (ELPD)  allow  a  higher  penetration  of  renewable  energies 
as  a  need  for  compliance  with  the  Kyoto  directives.  Moreover, 
their  lower  fuel  and  operational  costs  contribute  to  a  higher 
performance  in  terms  of  system’s  economy.  On  the  other  hand  in 
the  FPD  scenario,  higher  electricity  demand  levels  are  expected 
thus  imposing  the  requirement  for  intensive  use  and  higher 
capacity  orders  from  both  renewable  and  conventional  plant 
types.  Given  the  restricted  energetic  potentials  of  the  renewable 
sources,  the  conventional  plants  should  cover  higher  demand 
percentages  in  the  expense  of  the  aggregate  system’s  economic 
performance,  especially  due  to  their  higher  investment  and 
emission  costs.  In  both  case  studies  a  significant  improvement 
with  time  in  the  C02  emissions  per  MWh  produced  is  expected, 
compared  to  the  currently  recorded  values,  thus  indicating  an 
optimised  pattern  in  C02  emissions  and  therefore  to  the 
corresponding  emission  trading  costs. 

The  results  indicated  that  the  model  could  become  a  useful  tool 
for  energy  policy  advisors  and  authorities,  state  regulators  and 
decision  makers.  Potential  private  investors  and  electricity-related 
business  planners  could  also  benefit  from  the  tool  and  its  results 
since  they  may  have  an  overview  of  the  future  of  different 
investment  options  and  an  indication  of  the  shift  in  financial 
attractiveness  of  the  various  technologies  due  to  the  learning  curve 
effects.  As  long  as  they  may  engage  into  the  energy  market  arena 
they  could  identify  their  optimal  investing  options. 


Appendix  A.  Current  Greek  legal  status  notes 

In  Greece  the  Law  2773/99  constitutes  the  basic  legal  background. 
Two  companies,  the  Regulatory  Authority  of  Energy  (RAE)  as  well  as 
HTSO  S.A.  (or  the  so-called  “Hellenic  Transmission  System  Operator” 
or  even  simpler  “Operator”)  have  been  created.  HTSO  is  the  company 
that  handles  the  operation  of  the  Hellenic  Transmission  System  of 
Electric  Energy.  These  two  companies  are  the  basic  key  players  of  the 
deregulated  electricity  market.  The  main  key  players  of  the  Greek 
Electricity  market  are: 

-  RAE  is  an  independent  authority  that  manages,  suggests  and 
promotes  the  existence  of  equal  opportunities  and  fair  competi¬ 
tion  and  awards  the  operation  license  to  producers,  providers 
and  all  the  other  players,  related  to  the  market. 

-  HTSO  S.A.  is  a  company  with  a  twofold  role:  The  first  role  is  the 
one  being  played  by  P.P.C.  in  respect  to  the  Transmission 
System:  it  ensures  the  balancing  between  production  and 
consumption  and  that  the  electric  energy  is  provided  in  a 
reliable,  safe  and  in  terms  of  quality  acceptable  way  in  terms  of 
quality.  The  second  role  of  HTSO  is  to  settle  the  market,  in  other 
words  to  act  like  an  energy  stock  market  that  arranges  on  a 
daily  basis  who  owns  to  whom.  HTSO  does  not  provide  electric 
energy  and  all  the  basic  exchanging  relations  that  take  place 
are  bilateral  ones  between  producers/providers  and  their 
customers. 

-  The  Public  Power  Company  (P.P.C.)  is  one  among  a  number  of 
companies  that  will  be  created  in  the  field  of  electric  energy.  A 
stock  market  pictorial  description  of  the  roles  could  be,  P.P.C.  is 


an  admitted  company,  HTSO  is  the  stock  market  and  RAE  is  the 
Capital  Market  Committee. 


Appendix  B.  The  mathematical  background  for  the  forecasting 
tool 

Let  us  consider  the  following  general  stochastic  differential 
equation  (SDE): 

dXt  =  F(t,Xt)dt  +  G(t,Xt)dWt  (1*) 

Where: 

•  X  is  a  state  vector  of  process  variables  (for  example,  fuel  or 
electricity  prices)  to  simulate. 

•  W  is  a  Brownian  motion  vector  whose  differential  oscillates  in  a 
range  of  values  generated  by  normal  distribution. 

•  F  is  a  vector-valued  drift-rate  function. 

•  G  is  a  matrix-valued  diffusion-rate  function. 

The  drift  and  diffusion  rates,  F  and  G,  respectively,  are  general 
functions  of  a  real-valued  scalar  sample  time  t  and  state  vector  Xt.  The 
Eq.  (1*)  is  useful  in  implementing  derived  classes  that  impose 
additional  structure  on  the  drift  and  diffusion-rate  functions.  The 
derived  classes  of  models  used  in  this  study  are  the  following. 

GBM:  The  Geometric  Brownian  Motion  Class  whose  general 
equation  is: 

dXt  =  fjL(t)Xtdt  +  D(t,Xt)V(t)dWt  (2*) 

Where: 

•  Xt  is  a  state  vector  of  process  variables. 

•  /z(t)  is  the  mean  drift  function  derived  from  the  historical  data  as 
a  mean  difference. 

•  D(t,  Xt)  is  a  diagonal  matrix- valued  function.  Each  diagonal 
element  of  D  is  the  corresponding  element  of  the  state  vector 
raised  to  the  corresponding  element  of  an  exponent  (o'(t)),  which 
is  a  vector-valued  function:  D(t,X“(t)). 

•  V(t)  is  a  matrix- valued  volatility  rate  function. 

•  dW(t)  is  a  Brownian  motion  vector  (noise)  differential  which  is 
equal  to  e  Vdt  and  e  e  N( 0,1 ). 

MR:  The  Hull-White/Vasicek  (HWV)  Gaussian  Diffusion— Mean 
Reverting  Class: 

dXt=S(t)[L(t)-Xt]dt+V{t)dWt  (3*) 

Where: 

•  Xt  is  a  state  vector  of  process  variables. 

•  S(t)  is  the  reverting  speed  function  derived  from  the  historical 
data. 

•  L(t)  is  the  expected  long  run  mean  function. 

•  V(t)  is  a  matrix- valued  volatility  rate  function. 

•  dW(t)  is  a  Brownian  motion  vector  (noise)  differential  which  is 
equal  to  e  Vdt  and  e  e  N( 0,1 ). 
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